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EFFECT OF ENVIRONMENTAL TORQUES ON SHORT-TERM 
ATTITUDE PREDICTION FOR A ROLLING-WHEEL 
SPACECRAFT IN A SUN -SYNCHRONOUS ORBIT 

By Ward F. Hodge 
Langley Research Center 

SUMMARY 

The change in attitude produced by the action of environmental disturbance torques 
was investigated for a hexagonal cylinder spacecraft in a nearly polar 500-kilometer cir- 
cular orbit. A numerical evaluation and an analysis of the separate and combined effects 
of five such torques, arising from gravitational, magnetic, and surface pressure forces, 
were performed. The resulting perturbations were determined from time histories of 
the basic attitude state variables generated by numerical integration of the spacecraft 
equations of motion covering about one-third of the orbital period. Similar data were 
generated for examining linearity of the torque perturbations and sensitivities to system 
design changes and torque coefficient variations. 

The torque perturbations proved to be very small and exhibited linearity such that 
linearized equations of motion yielded accurate results over short periods and the sepa- 
rate perturbations contributed by the individual torques were additive in the sense of 
superposition. Linearity of the torque perturbations was not affected by moderate system 
design changes and persisted for torque-to-angular momentum ratios up to 100 times the 
expected 10 nominal value. As these conditions include many possible applications, 
similar linear behavior might be anticipated for other rolling-wheel spacecraft. 

The principal perturbations consisted of a relatively small spin decay of about 1° 
caused primarily by eddy-current torque, and a slight spin-axis drift of some 70 arc- 
seconds due mainly to solar radiation pressure and magnetic -moment torques. Because 
of the small 10 “7 torque-to-angular momentum ratio, the affected torque terms in the 
spacecraft equations of motion remain so small that Euler -angle singularities are unlikely 
to be approached closely enough to cause difficulties. In regard to attitude prediction, the 
environmental torques were estimated to have a much lesser effect than the propagation of 
state estimation residuals. 


INTRODUCTION 


The use of spin- stabilized satellites for conducting measurements research, in 
fields such as earth resources and meteorology, generally requires precise spacecraft 
attitude information over extended periods. For such applications, the line of sight of an 
instrument axis fixed in the spacecraft must be known within about 15 arc-seconds during 
data-gathering periods to obtain results having the desired spatial resolution. This prob- 
lem was investigated as part of recent studies of spacecraft configurations and instru- 
mentation for high resolution measurements of the earth's infrared horizon profile con- 
ducted or sponsored by the NASA Langley Research Center. 

The stringent accuracy required for this type of mission necessitates frequent 
updating of the spacecraft attitude to prevent the propagation of state estimation residuals 
and dynamical modeling uncertainties from causing the allowable pointing error to be 
exceeded. Attitude prediction or state transition over long periods is therefore not fea- 
sible, and a short-term approach must be implemented. The scheme employed for the 
horizon measurements studies assumes the necessary updating can be satisfactorily 
accomplished by means of precise attitude estimation. Based on operational limitations 
which affect the availability of attitude determination data, the maximum interval between 
updates for which attitude prediction would be required is about one-third of the space- 
craft orbital period. 

In seeking a suitable analytical procedure for achieving the required attitude accu- 
racy and handling problems such as Euler-angle singularities, two alternate choices for 
the basic attitude state variables were studied independently and concurrently. Because 
of their influence on attitude prediction, an important task in each case was to model and 
evaluate with care the effects of significant disturbance torques. The systems study of 
reference 1 and the excerpt from it given in reference 2 contain torque calculations based 
on the use of state variables comprised of the body-axis components of the spacecraft 
angular velocity, and aircraft -type Euler angles analogous to yaw, pitch, and roll. While 
the computations from both studies were essentially in agreement, a more extensive eval- 
uation of the torques and their effects was performed by using the angular momentum 
formulation described in the present paper. The results and conclusions obtained are 
more general than those of references 1 and 2, and the analysis presented includes lin- 
earity of the torque perturbations, sensitivities to system design changes, and Euler-angle 
singularities. In summarizing the overall investigation, the two approaches were com- 
pared and the use of time-averaged equations of motion such as those of reference 3 were 
discussed. 
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SYMBOLS 


A area of spacecraft panel exposed to surface pressure, centimeters^ 

Ax,Ay,A 2 aerodynamic pressure torque coefficients, dyne-centimeters 
a altitude, kilometers 

B geomagnetic field intensity or flux density, teslas 

b height of spacecraft side panels, centimeters 

C(x 0 ,<i> 0 ,£ 0 ) spacecraft panel orientation matrix (eq. (B3)) 
c spacecraft eddy-current coefficient, meters4/ohm 

E (<p,6,\p) spacecraft Euler-angle matrix (eq. (3)) 
e eccentricity of spacecraft orbit 

ej elements of E (<p,6,if/) 

F surface pressure force, dynes 

G(of,0) geomagnetic field matrix (eq. (A5)) 

gj^jh^ 11 Schmidt-normalized Gaussian coefficients, nanoteslas 
H(£,t) spacecraft angular momentum vector orientation matrix (eq. (2)) 

h spacecraft angular momentum, dyne-centimeter- seconds 

I^lyjlj, spacecraft principal moments of inertia, gram -centimeter s^ 
i inclination of spacecraft orbit, degrees 

L geocentric true longitude of sun on ecliptic, degrees 

l length of spacecraft solar cell panels, centimeters 


3 



M 


spacecraft residual magnetic moment, ampere(turn) -meters^ 
iiq spacecraft mean orbital motion, radians/second 

P^ n (0) Schmidt-normalized associated Legendre function 

p unit vector along incident direction of a surface pressure 

q Q aerodynamic pressure, dynes/centimeter2 

R^RyjRg solar radiation pressure torque coefficients, dyne -centimeters 
r radius, centimeters 

r a mean radius of earth, kilometers 

r Q geocentric distance of spacecraft, kilometers 

r^ equatorial radius of earth, kilometers 

S intensity of solar radiation pressure, dynes/centimeter^ 

s fraction of p reflected specularly (eqs. (B8)) 

T environmental torque, dyne- centimeters 

t time, seconds 

u argument of latitude of spacecraft orbit, degrees 

v 0 linear velocity along spacecraft orbital path, kilometers/second 

W(f2,i,u) orbit orientation matrix (eq. (1)) 
w width of spacecraft panels, centimeters 

x,y,z rectangular coordinates and spacecraft principal axes 

a geocentric right ascension, degrees 
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effective spacecraft angle of attack, degrees 
solar elevation phase angle, degrees 

direction cosines of r" 0 relative to spacecraft principal axes 

geocentric declination, degrees 

obliquity of ecliptic, degrees 

surface pressure angle of incidence, degrees 

geocentric colatitude, degrees 

geocentric longitude, measured east from Greenwich meridian, degrees 
spacecraft panel orientation angles, degrees 
universal gravitational constant, kilometer s^/second^ 
surface pressure azimuth relative to a spacecraft panel, degrees 
geocentric right ascension of h, degrees 

fraction of incident solar radiation reflected by A (eqs. (B8)) 
atmospheric density, kilogram s/meter ^ 

aerodynamic pressure surface reflection coefficients for normal and 
tangential momentum exchange 

geocentric codeclination of h, degrees 

sum of Euler angles <f> and i//, degrees 

spacecraft Euler angles, degrees 

geocentric right ascension of ascending node of spacecraft orbit, degrees 
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spacecraft angular velocity, radians/second 


U) 

a) axial rotation rate of earth, radians/second 

© 

Subscripts: 

G Greenwich meridian 

h angular momentum 

i inertial 

j,k matrix element indices 

N normal component 

n degree of spherical harmonic expansion 

o orbital, initial, or reference value 

T tangential component 

x,y,z spacecraft principal axis components 

0,A,r geocentric spherical polar components 

Superscript: 

m order of spherical harmonic expansion 

Notation: 

f transformed, deviated, or compared quantity 

time derivative 
~ unit vector 

— vector 
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centroidal distance 


•= centroidal vector 

-1 matrix inverse 

A residual or differential value 

sec arc-second 

ANALYTICAL FORMULATION 

The angular momentum formulation used in the study is described in two steps 
which respectively define the spacecraft system and the dynamics of the rigid-body 
motion. The first concerns the spacecraft configuration, orbit, and attitude, and the sec- 
ond, the equations of motion and modeling of the environmental disturbance torques. 

Spacecraft Configuration 

The basic spacecraft assumed for both studies was a rolling-wheel vehicle in the 
form of a right hexagonal cylinder having six solar cell panels attached to one base. As 
depicted in figure 1, this configuration consisted of two hexagonal end panels of width 2w 
across corners, six side panels of width w and height b, and six solar cell panels of 
negligible thickness, width w, and length l. The spacecraft principal axes x, y, and 
z and center of mass were assumed to have the orientations indicated in the sketch. The 
dimensions w, b, and l and other pertinent spacecraft data are listed in table I. 



Figure 1.- Hexagonal spacecraft configuration. 
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Spacecraft Orbit 


A geocentric equatorial coordinate system x^y^zj, having Xj toward the vernal 
equinox and z^ along the earth's spin axis as shown in figure 2, was used as the inertial 



Figure 2.- Spacecraft orbit. 


reference for both the spacecraft orbital and rigid-body motions. The orientation of the 
spacecraft orbit relative to these axes is defined by the three angles £2, i, and u as 
indicated. Let x 0 , y 0 , and z 0 be a set of rotating axes alined with the orbit plane such 
that z 0 coincides with the normal to the orbit plane and x 0 is along the local vertical. 
The transformation between the orbital and inertial coordinate systems is then 


l x o) 


h\ 

y 0 

= W(£2,i,u) 

?i 


f 

W 


where 


W(0,i,u) 


cos u cos £2 - cos i sin u sin £2 
-sin u cos £2 - cos i cos u sin £2 
^ sin i sin £2 


cos u sin £2 + cos i sin u cos £2 
-sin u sin £2 + cos i cos u cos £2 
-sin i cos £2 


sin u sin i\ 
cos u sin i 
cos i 


The relevant physical data and values assumed for the orbital elements may be found in 
tables II and m, respectively. 
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Spacecraft Attitude and Pointing Direction 

The manner in which the angular momentum state variables define the spacecraft 
attitude is described. Then the angular coordinates of a spacecraft axis, chosen to repre- 
sent the line-of-sight or pointing-direction reference for onboard instrumentation, are 
expressed as functions of these variables. 

State variables . - The angular momentum state variables express the spacecraft 
attitude in terms of five angles and the magnitude of the angular momentum vector. The 
first two angles £ and t define the inertial direction of the angular momentum vec- 
tor K, and the remaining three angles <p, 6, and if/ give the orientation of the space- 

craft principal axes x, y, and z relative to h. Figure 3 illustrates the definitions of 



Figure 3-- Orientation of spacecraft angular momentum vector. 


£ and r, and an angular momentum coordinate system xj^yj^Zh chosen such that z ^ 
coincides with K. The transformation of these axes to inertial reference is 


l x h\ 


! x i\ 

yh 

= H(£,t) 

n 

Wi 


W 


where 


H(|,r) = 


cos £ 

-sin £ cos t 
^ sin £ sin t 


sin £ 

0 

cos £ COS T 

sin t 

-cos £ sin r 

COS T 


( 2 ) 
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The Euler angles <£, 8, and x j/, which relate the principal axes to the system 

x h>yh> z h in ^e manner shown in figure 4, complete the definition of the spacecraft 


z h 



attitude in terms of the angular momentum state variables. The x-, y-, and z-axes and 
x h -, Yh - ? and z^-axes are related by 




/ x h\ 

y 

= E 

yh 

W 


w 


where 



j cos i p cos $ - cos 6 sin (p sin i fy 

cos 0 sin 0 + cos 0 cos 0 sin 0 

sin 0 sin $\ 

E(<M,«/') = 

-sin cos 0 - cos 0 sin 0 cos 0 

-sin 0 sin 0 + cos 0 cos 0 cos 0 

cos 0 sin 9 


^ sin 6 sin 0 

-sin 9 cos 0 

cos 9 


Pointing direct ion. - The geocentric right ascension a and declination 6 of the 
principal axis x were used to define the spacecraft pointing direction as indicated in 
figure 5. As the instrument line of sight was intended to be perpendicular to the spin 
axis z, this choice provides a suitable basis for determining the net change in pointing 
direction caused by the torques. Expressions for evaluating a and 6 in terms of the 
state variables t, cp, 8, and i// were obtained by equating the direction cosines of 
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Earth spin axis 




Figure 5.- Spacecraft pointing direction. 

the x-axis as defined by a and 6 in figure 5 with those resulting from the total trans- 
formation given by equations (2) and (3). The desired relationship is 

eu cos £ - sin ^e^2 cos 7 ~ e 13 s i n T ) 
eu sin £ + cos cos 7 ~ e 13 s i n T ) 

ej2 sin r + ei3 cos r 

from which 


(4) 


where e-Q, e-^, an< 4 e i3 and the right-hand members of equations (4) are the ele- 
ments of the first rows of E and the product EH, respectively. 

Equations of Motion 

Although the customary form of Euler's dynamical equations provides a suitable 
description of the spacecraft rigid-body motion, equivalent first-order expressions for 
the rates of the six state variables were used instead. These equations are somewhat 
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e ll s * n £ + cos £ ( e 12 cos T " e 13 s * n T ) 
en cos | - sin £ (e^2 cos 7 ~ e 13 si n T ) 

sin 6 = e^2 sin 7 + e i3 cos r 




more convenient for numerical integration and differ only slightly from those of refer- 
ence 4 which extend earlier formulations (for example, see ref. 5) to include torques. 
The expressions for r, h, <£, 9, and 4/ used in the study were obtained as fol- 

lows. By starting with the equation for the torque about the center of mass of a rigid 
body 


T = h + (uxh) 
or 


l T *\ 




^ ^yhz ” w zhy ^ 

T y 

= 

hy 

+ 

-u> x h z + w z h x 

1 



\ KZ / 


\ w x hy - coyh X J 


j 


the principal axis components of h 


h x = I x u> x = h s in ^ sin ^ 
hy = IyCi>y = h sin 9 cos 4 / 
h z = I z cu z = h cos 9 


(6) 


may be used to eliminate the components of 3 from equations (5). Substitution of equa- 
tions (6) and their time derivatives into equations (5) gives 

T x = h sin 9 sin 4 / + h(fi cos 6 sin ip + 4 1 sin 0 cos ’Z') + h2 sin 9 cos 9 cos 4 / - j-j 

Ty = h sin 9 cos 4 / + h((9 cos 9 cos \j/ - sin 9 sin xp) - h^ sin 9 cos 9 sin \ p ~ j~j 

T z = h cos 9 - h 9 sin 9 + h2 sin^ 9 sin j// cos 4/ \~ - -r- ) 

V x l V 

from which 


h = (t x sin 4/ + Ty cos sin 9 + T z cos 9 


9 = h sin 9 sin 4/ cos 4/ - — ) + 


r x hi 


^T x sin 4* + Ty cos \p) cos 

i = h cos 6 (f - 5isi Jfc . SSpk\ + ;-J— (T x cos * - T y sin 
\I Z I x Iy J h sin 9\ A y 


9 - T z sin 9 

*) 
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The remaining three equations were obtained in a similar manner by substituting 
another definition of the components of a> into equations (5). Reference to figure 4 
shows that 


cf> sin 6 sin if/ + 6 cos iM 


Wy = E {<p,6,\p) | sin r + $ sin 6 cos - 6 sin 

10 O z COS T \(p COS 9 + 4/ j 


where the first term on the right-hand side accounts for the angular velocity of h. By 
substituting equations (7) into equations (5) and performing the indicated algebraic 
reductions 

T x = h sin 9 sin ^ + h^-ej^r + e lll sin T ) 

T y = h sin 9 cos ip - h^e22^ - sin T ) 

T z = h cos 9 + h sin 9 (r cos (p + £ sin t sin <pj 
which may be solved for 

t = - - — f — -h cos 9 cos <p + sin 9 sin <p (T v cos il/ 
h sin 0l_ V x 

and 

£ = h S i n \ s ' fn cos ® sin ^ + sin 0 cos 0 ( Tx cos ^ ” T y sin + T Z sin 

where h is given by the expression previously obtained. Substituting oj z from equa- 
tions (7) into the third member of equations (6) and solving for <f> gives 

4> = - cQ ^ - g - £ cos r - tan 6 (t sin (f> - £ sin r cos <pj 
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Thus the complete set of equations is 


h = (t x sin i// + Ty cos i j/j sin 6 + T z cos 6 

4 = h s ^ n T jcos <P (T X cos ^ - T y sin i - cos 6 sin <p ^T x sin 4/ + T y cos i//) + T z sin 0 sin 

f = -ijsin <p (t x cos if/ - T y sin + cos 6 cos (j> (t x sin 4/ + T y cos i//) - T z sin 0 cos <^J 

4 = h + SS^j . £^(t x cos 4, - Ty sin *) - 4 COS r 

0 = h sin 0 sin 4^ cos ^ (t x sin 4 / + T y cos xf) cos 0 - T z sin 0 J 


> (8) 


Environmental Torques 

Of the possible sources of disturbance torques, only those arising from gravita- 
tional, magnetic, and surface pressure forces appeared large enough to affect the space- 
craft attitude. The resulting analytical model consisted of five environmental torques 
due to gravity gradient, magnetic moment, eddy current, solar radiation pressure, and 
aerodynamic pressure. All other torques were estimated to be either negligible over the 
short intervals being considered or rendered insignificant by a systems operation con- 
straint that all energy radiators, dampers, and control systems be deactivated during 
data-gathering periods for which precise attitude information is required. Torque 
sources such as mass expulsion and momentum transfer were further minimized by using 
magnetic coils in place of gas-operated and momentum exchange control devices. 

For the spacecraft inertia configuration assumed in table I, the spherical earth 
gravity- gradient torque model given in reference 6 was considered adequate. The 
torques due to interaction with the geomagnetic field were evaluated from the well-known 
equations appearing in reference 7 and the field computation procedure described in 
appendix A. As the torques arising from surface pressure forces are highly dependent 
on the geometry and orientation of the spacecraft surfaces, the solar radiation and aero- 
dynamic pressure torque equations usually must be formulated to suit the application. 

The equations developed in appendix B for the hexagonal configuration are based on the 
surface pressure force models contained in references 8 and 9. 
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Gravity gradient. - According to the spherical earth model, the principal axis com- 
ponents of the gravity-gradient torque are 


Tx = 3n 0 ^I z - I y )y y y z 
Ty = 3n 0 ^I x - I z )y x y z > 
Tz = 3n 0 ^^Iy - I x ) y x y y 


( 9 ) 


where n 0 = i/p/r 0 3 is the spacecraft mean orbital motion (see table ID) and 


/>x' 


M\ 

r y 

= EHW-l 

0 

W 

1 

\°) 


are the direction cosines of the gravitational attraction force with respect to the x, y, 
and z principal axes as obtained from equations (1) to (3). 

Magnetic moment .- The torque caused by the interaction of the residual magnetic 
moment of the spacecraft with the geomagnetic field is 


T = M X B 
or 

T x = MyB z - M z By l 
Ty = Mg-Bjj - M X B Z 
Tz = Mx®y “ MyBjj^ 


(ID 


where the values assumed for M x , My, and M z are listed in table I, and Bjj, By, 
and B z may be evaluated by means of the procedure described in appendix A. 
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Eddy curre nt.- Assuming the spacecraft hull to be approximated as a thin spherical 
shell, the torque due to eddy-current circulation is 

T = c(w x b) X B 
or 

T x = -cjo> x ^By2 + B z ^ - B^ByWy + B z o> z jj 

Ty = -c|a;y^B x ^ + B z ^) - By^o^ + B z u> z ^ 

T z = -c |^ , z(b x ^ + By^ - + ByCOyjj 

where the value used for the eddy-current coefficient c is also listed in table I and co x , 
coy, and co z may be computed from either equations (6) or (7). The adequacy of equa- 
tions (12) for the hexagonal configuration was demonstrated in references 1 and 2 where 
a more elaborate cylindrical shell model was shown to be unnecessary. 

Solar radiation pressur e. - Based on the evaluation performed in appendix B, the 
principal axis components of the solar radiation pressure torque were assumed to be 
reasonably approximated as 


T x = -Rx sin ' 0) 
T y = -R y cos (# - p) l 

T z = R z sin 6 (# - [3 ) 


(13) 


where $ = <p + \p and j3 is a solar elevation phase angle (see appendix B). The values 
derived for the torque coefficients R x , R y , and R z may be found in table I. 

Aerodynami c pressure.- A similar evaluation (see appendix B) yielded 


T x = A x | sin Q! 0 j cos u cos (<3> - u) 

Ty = -Ay | sin cv 0 1 cos u sin ($> - u) \ 
T z = -A z sin a Q cos u sin 6 (<3? - u) 


(14) 


for the principal axis components of the aerodynamic pressure torque where a 0 is the 
effective angle of attack and A x , Ay, and A z are also derived torque coefficients (see 
table I). 
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RESULTS AND ANALYSIS 


A numerical evaluation and an analysis of the torqued spacecraft motion, based on 
equations (8) to (14), were performed for the nominal conditions and data listed in tables I 
to III. The procedures employed for numerically integrating these equations and deter- 
mining the total instrument-axis pointing error are described; then the deviations or 
perturbations caused by the environmental torques are analyzed. These quantities were 
calculated from time histories of the state variables covering 2000 seconds or about a 
third of the orbital period, which was the maximum prediction interval required as stated 
in the introduction. Typical results for only one 2000-second period are presented since 
time histories obtained for nominal and worst case conditions indicated that the nature of 
the torque perturbations should remain essentially the same for all such intervals. 

Lastly, the angular velocity and angular momentum formulations were compared, and 
the use of time-averaged equations of motion was discussed. 

Approach 

The computational considerations addressed in generating the necessary time 
histories were the choice of numerical method, comparison of direct and perturbation 
solutions of equations (8), effect of singularities on the use of angular momentum state 
variables, and calculation of the instrument-axis pointing error. 

Numerical integration .- As equations (8) are of first order and the torques are very 
small, the fourth-order Runge-Kutta scheme was chosen as the integration method. The 
integration step was determined by progressively reducing the step size until successive 
solutions differed by less than the equivalent of 2 arc -seconds in all six state variables 
over the full 2000-second integration period. This level was reached when the step size 
was decreased from 0.1 to 0.01 second. A step size of 0.1 second thus appeared adequate 
and was chosen in anticipation that the total integration error due to truncation and round- 
off probably would not exceed a few arc-seconds. 

In order to determine whether a perturbation approach would be advantageous, time 
histories generated from equations linearized about the untorqued solution of equations (8) 
were compared with corresponding ones obtained by direct integration. Both methods 
were found to give essentially the same results over the entire 2000 -second period with 
no rectification of the linearized solution, thus indicating that the small torque perturba- 
tions were accurately approximated by using the linearized equations. 

Euler -angle singularities . - Inspection of equations (8) shows that singularities will 
occur if either 6 or r is zero. The latter was of no concern as attitude control will 
constrain t to within 5° of the nominal 97.72° value listed in table in, but the 
^-singularity could occur since the spacecraft was intended to spin without coning. How- 
ever, no difficulty was encountered for the 0.66° coning angle assumed for the study. 
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This result was verified by repeating the computations using the angular velocity state 
variables of references 1 and 2 which avoid all singularities for the conditions being 
considered. The two sets of values for the instrument -axis coordinates a and 6 
were identical within a few hundredths of an arc-second. 

The reason small values of 6 caused no trouble stems from the smallness of the 
ratio of T x cos if/ - Ty sin if/ to h in the fourth and sixth members of equations (8). 
Comparison of this ratio, which was the order of 10“^ for the data given in tables I and 
HI, with sin 6 shows that 6 could approach zero within a few arc-seconds before the 
affected torque terms become even as large as the untorqued ones. The likelihood of 6 
becoming this small appears remote because the minimum initial coning angle achievable 
with existing nutation dampers is limited by factors such as tuning errors and residual 
spacecraft imbalance. Further inspection of equations (8) indicates that t also could 
be within a few arc-seconds of zero without causing difficulty for the same reason. Thus, 
the use of angular momentum state variables for the horizon measurement spacecraft 
does not appear to be restricted by either singularity and need not be avoided. 

Instrument -axis pointing error .- The uncertainty in the instrument line of sight, 
which ultimately must be determined, was considered to be suitably represented by the 
total change in the pointing direction of the assumed instrument axis. This quantity gives 
a direct measure of the total line-of-sight error caused by the torques and clearly indi- 
cates the relative significance of each torque. The pointing-direction error was defined 
to be the total change in the inertial orientation of the assumed instrument axis (see 
fig. 5). This angle is denoted as Ax in figure 6 and may be readily calculated from 



Figure 6.- Total pointing-direction deviation. 
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equations (4) and either the law of cosines formula 


cos Ax = sin 6 sin 6' + cos 6 cos 6' cos (a - a') 
or the approximation 

Ax ® \/(Aq! cos 6)2 + (A6)2 


(15a) 


(15b) 


where 


A a = a - a' 


A6 = 6 - 6' 


and (a, 5) and (a', 6') refer to the original and deviated positions of the instrument 
axis, respectively. 


Computational Results 

Six torque cases requiring seven sets of time histories were used to determine the 
separate and combined effects of the five torques. The results are presented as plots of 
the pointing direction and state variable deviations obtained as the differences between 
appropriate pairs of time histories. Similar data were generated for investigating linear 
characteristics of the torque perturbations, and sensitivities to torque coefficient varia- 
tions and system design changes. Discussion of the results is concluded with an evalua- 
tion of the effects of the torques on attitude prediction. 

Time-history plots .- The changes in the pointing direction and state variables 
obtained for each of the six torque cases are plotted as figure 7 and the corresponding 
torques as figure 8. Parts (a) to (f) of each figure present the six sets of results in the 
following order: combined torques, gravity gradient, magnetic moment, eddy current, 

solar radiation pressure, and aerodynamic pressure. The total deviations produced by 
the simultaneous action of all five torques, plotted as figure 7(a), represent the differ- 
ences between an untorqued solution of equations (8) and one having the five torques 
included. The remaining parts of figure 7 present similar plots of the individual effects 
of each torque, and the pointing direction error Ax for all six parts was computed from 
equations (4) and (15). The three frequencies evident in both figures correspond to those 
of 0, 0, and u, which had respective periods of about 20, 128, and 5680 seconds for the 
conditions assumed in tables I to m. Quantities which were zero account for the curves 
missing from some parts of figures 7 and 8. 
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(f) Aerodynamic pressure. 


Figure 7*- Concluded. 
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(c) Magnetic moment. 
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(d) Eddy current. 


Figure 8.- Continued. 










The results presented as figure 7 indicate that the principal effect of the torques 
was to despin the spacecraft in nearly an exponential manner with little change in the 
inertial direction of the principal spin axis. Reference to the plotted curves shows that 
the instrument-axis pointing direction was progressively altered by roughly 1° over 
2000 seconds or about 1/36 000 of its untorqued spin motion. Figure 7 further indicates 
that Ax was mainly the result of a reduction in spin rate, and was approximately equal 
to the algebraic sum of A (f> and Ai {/. This reduction is also evident from the curves 
for Ah which exhibit a proportional decrease in angular momentum. Because the small 
coning angle of 0.66° assumed in table III was almost unchanged by the torques (see 
fig. 7), the principal spin axis remained nearly alined with IT so that the slight spin 
axis drift of some 70 arc-seconds is exhibited mainly by the curves for A| and At. 
These results show the spin decay to be primarily due to eddy-current torque, and the 
slight spin axis drift, to solar radiation pressure and magnetic-moment torques. The 
perturbations caused by environmental torques thus appeared to be very similar to those 
expected for most rolling-wheel spacecraft. 

Linearity . - The linear nature of the torque disturbances mentioned previously, and 
referred to as an additive property in references 1 and 2, was also evident in the compu- 
tations being discussed. Addition of the separate perturbations contributed by each torque 
in the sense of superposition yielded very nearly the same results as those obtained by 
considering all five simultaneously. The smallness of the torques is clearly the reason 
linearity is exhibited by the perturbations . 

In order to determine the extent or region of linearity, additional solutions to equa- 
tions (8) were generated with the combined total torques multiplied by factors of 2, 10, 

100, and 1000. These factors represent torque levels having respective potentials for 
angular momentum change ranging from 10“^ to 10“^ of the nominal magnitude for h 
given in table EH. The corresponding amounts by which superposition failed to hold were 
respectively 0.05, 0.25, and 25 arc-seconds for the first three factors and about 13° for 
the fourth, thereby indicating fairly good agreement for torque amplitudes up to 100 times 
their nominal expected levels. The linear region thus defined should also include many 
other spin-stabilized spacecraft as the torque-to-angular momentum ratio T/h probably 
will not exceed this 10"® linearity limit in most cases. For example, rough calculations 
show this ratio to be only 10“6 to 10“^ for spacecraft such as Tiros, Tiros M, DME-A, 
and even Telstar and Syncom. 

Sensitivity . - The question of sensitivity to possible variations or errors in torque 
coefficients was easily resolved as linearity also causes the perturbations to be in direct 
proportion to the torques that produce them. A 10-percent change in the eddy-current 
coefficient, for example, would be expected to alter the amplitudes of the curves plotted 
as figure 7(d) by about 10 percent. This property is of interest for estimating the effects 
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on attitude prediction of variations in some torque coefficients that may occur between 
the dark and sunlit regions of the spacecraft orbit. 

In regard to system design factors, configuration changes in the spacecraft axial- 
to -transverse inertia ratio and asymmetry were of concern because of their possible 
effect on the dominant eddy-current torque. For inertia ratios up to 1.5 and asymmetries 
as large as 10 percent, which includes most rolling-wheel configurations of interest, the 
torque perturbations exhibited the same linear properties and were very similar to those 
presented as figure 7. Based on the evaluations performed in reference 1 and appendix B, 
geometrical factors such as spacecraft panel arrangement should have little effect on the 
eddy-current and surface -pressure torques unless fairly extensive changes are made or a 
paddle-wheel configuration is used. Halving and doubling the spin rate also had minor 
effects on the results, which was as expected since T/h for the eddy-current torque 
does not change with the magnitude of c o (see eqs. (6), (8), and (12)). As neither the 
overall nature or linearity of the torque perturbations were altered by moderate system 
design changes, similar results might be anticipated for other rolling-wheel spacecraft. 

Attitude prediction. - As mentioned in the introduction, attitude prediction accuracy 
was a primary concern in investigating the effects of environmental torques. The pre- 
diction error incurred by ignoring a given torque, which can be readily estimated from 
figure 7 or similar results, was used in references 1 and 2 as a criterion for determining 
whether any previously modeled torques could be omitted from the computations. This 
criterion serves to verify which torques must be included but does not account for the 
effects of possible uncertainties in their coefficients. Assuming the maximum such 
uncertainty to be 10 percent, reference to figure 7 indicates that only the eddy-current 
coefficient is likely to cause the specified error bound to be exceeded. 

Because the torques were very small, a comparison of their effects on prediction 
accuracy with those arising from the propagation of state estimation residuals was of 
interest. The achievable attitude determination accuracy for the horizon measurements 
spacecraft was estimated to be about 10 arc-seconds, for which equations (8) would be 
roughly 100 times more sensitive to state variable residuals than to all five torques com- 
bined. The resulting attitude error would not be 100 times larger as the effects of some 
residuals tend to cancel, but still much greater than that incurred by ignoring the torques 
completely. While state estimation accuracy thus appears to be the more critical problem 
and places a greater restriction on attitude prediction, the torques are not negligible and 
should be included. 


Comparison of State Variable Systems 

The overall results of both studies indicated little reason to favor either choice of 
state variables for the type of application considered. Although the use of somewhat 
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different conditions and alternate surface pressure torque models prevented exact com- 
parisons, the two sets of results were essentially in good agreement as mentioned in the 
introduction. The numerical differences stem mainly from the fact that the computations 
presented in references 1 and 2 are for worst case conditions, while figures 7 and 8 were 
generated for nominal ones. The exact comparisons made in conjunction with evaluating 
Euler -angle singularities further showed that both formulations yield almost identical 
results under the same conditions. A somewhat more compact analytical format is 
afforded by the angular velocity formulation, but the angular momentum state variables 
are simpler functions of time which permit better insight and physical interpretation of 
the rigid-body motion. With regard to time -averaging the equations of motion as dis- 
cussed subsequently, the angular momentum variables are much more convenient as the 
angular velocity equations of motion are awkward and impractical to average. 

An additional scheme of interest for evaluating the torque perturbations is the use 
of time-averaged equations of motion such as those of reference 3. While these equa- 
tions may cause some loss of accuracy, they simplify the attitude estimation process by 
greatly reducing the computational burden associated with their numerical integration. 

The application of this technique to equations (8) to (14) essentially removed the variations 
due to <p and \p such that the results nearly duplicated figure 7 as viewed with these 
harmonics deleted. Examination of figure 7 shows that the error introduced by time- 
averaging also differs for each torque and respectively ranges from about 3 to more than 
15 arc-seconds for the eddy-current and magnetic-moment torques. Depending on which 
torques are included in the computations, time-averaged equations of motion can thus be 
nearly an order of magnitude less accurate than indicated in references 1 to 3. 

CONCLUSIONS 

The results obtained from a numerical evaluation and an analysis of the effects of 
environmental disturbance torques on the attitude of a hexagonal rolling-wheel spacecraft 
indicate the following conclusions: 

1. The torque perturbations are very small and exhibit linearity such that linearized 
equations of motion yield accurate results over short periods and the separate perturba- 
tions contributed by each torque are additive in the sense of superposition. 

2. Linearity of the perturbations is not affected by moderate system design changes 
and should persist for torque -to -angular momentum ratios up to 100 times the nominal 
expected value, thus indicating that similar behavior might be anticipated for other 
rolling-wheel spacecraft as many possible applications are included by these conditions. 
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3. The principal perturbations in most cases will probably be a relatively small 
spin decay due to eddy-current torque and a slight spin-axis drift caused by solar radia- 
tion pressure and magnetic-moment torques. 

4. Because of the small torque-to-angular momentum ratio, Euler-angle singulari- 
ties should not cause computational difficulties or restrict the use of angular momentum 
state variables. 

5. The environmental torques will affect attitude prediction to a much lesser extent 
than the propagation of attitude estimation residuals in most cases. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 15, 1971. 
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APPENDIX A 


GEOMAGNETIC FIELD COMPUTATIONS 

In order to minimize errors in evaluating the magnetic torques, the magnetic field 
of the earth was computed from the spherical harmonic model and Gaussian coefficients 
g™ and h™ of reference 10 which includes terms through 10th order and degree. As 
the atmosphere is essentially nonmagnetic, the equations for the geocentric field intensity 
components given in reference 10 may be expressed in terms of the magnetic flux density 
as 


oo n 


B© = - Y, ^ (^j (Sn 1 cos mX + h n 1 sin mA 
n=l m=0 


vdP“(0) 

7 d© 


00 n n+ 2 
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=1 m=0 


°° n n+2 

B r = ^ ^ (n + 1 )(~j cos mX + h“ sin mAjP™(0) 
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where the Schmidt-normalized associated Legendre functions P™(©) and their deriva- 
tives may be evaluated from the recursive formulas 
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APPENDIX A - Continued 


which are much more efficient for digital computer programing than the familiar power 
series expansions for these functions. 

As equations (Al) give B©, B^, and B r as functions of r, 0, and X, the initial 

step in computing the desired principal axis components of the geomagnetic field was to 
express these coordinates in terms of the spacecraft orbital position. The necessary 
angular relationships for this purpose are illustrated in figure 9, which also indicates the 



Figure 9*- Geomagnetic field at spacecraft. 

spatial directions of B 0 , and B r . In terms of the notation of figures 2 and 9, 
expressions for evaluating 0 and X were obtained by means of the well-known 
formulas 


o; = £1 + tan"-*- (cos i tan u) 
0 = cos - -*- (sin i sin u) 


which appear in standard texts on celestial mechanics. With the geocentric right ascen- 
sion of the field point at the spacecraft determined from the first of equations (A3), the 
corresponding east longitude from the Greenwich meridian is 

X = a - otQ (A4) 

where 

a G = + a G 0 
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APPENDIX A - Concluded 


and Q!^ denotes the right ascension of the Greenwich meridian at a chosen epoch or 
initial reference time. 

After solving equations (Al) to (A4) for the proper values of B@, B^, and B r , the 

principal axis components B^, By, and B z may be obtained by means of the trans- 
formation 



= E m,r) G(a,Q) 



(A5) 


where E((p,d,>p) and H(£,r) are respectively defined in equations (3) and (2) and 


sin 0 cos a -sin a 
G(q!, 0) = sin 0 sin a cos a 
\cos © 0 


cos © cos a 
cos © sin a 
-sin 0 


gives the intermediate transformation of B@, B^, and B r to inertial axes. 
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APPENDIX B 


SURFACE -PRESSURE TORQUES 

The equations used for calculating the two surface -pres sure torques considered in 
the study were based on the solar radiation and aerodynamic pressure force models of 
references 8 and 9, respectively. These models express each force in terms of normal 
and tangential components and were chosen for convenience in obtaining general equations 
for any surface-pressure torque. The resulting expressions were applied to each of the 
two torques, then evaluated for the hexagonal spacecraft configuration. Because the 
resultant forces on the spacecraft panels act at the centers of pressure of the exposed 
surface areas, a partly computerized semigraphical procedure for locating the corre- 
sponding centroids was employed. This approach avoided time-consuming repetitive 
evaluations of surface integrals by taking advantage of the regular geometry of the panels 
and overall uniformity of the spacecraft rotational motion and led to efficient approximate 
expressions for these torques. 


General Formulation 

With reference to figure 10 let A be that area of a spacecraft panel exposed to a 
surface pressure. The resultant force F then acts through the center of pressure or 


x’ 
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APPENDIX B - Continued 


centroid of A as indicated. Relative to the center of the panel, which is the origin of 
the x'-, y'-, and z' -axes oriented such that y’ and z' lie in the plane of the panel with 
x' normal to it, the centroid of A is located by the vector 



where x' is zero in all cases. The incident direction of the surface pressure is 
coplanar with F and the x'-axis as indicated in figure 10 and is given by the unit vec- 
tor p the direction cosines of which are 


P = 


cos rj 

sin 77 cos v 
sin r\ sin v 


The resolution of F into components parallel to the x'-, y'-, and z’-axes is 


(Bl) 


F = 



where 


(B2) 


F x' ~ “ F N 
Fy, = -F t cos v 

F z t = -F-p sin v 

in which f N and F-p are, respectively, the normal and tangential components of F. 

The orientation of the panel relative to the spacecraft center of mass is defined by 
the distance 
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between the origins of x', y', and z' and the principal axes x, y, and z, and the 
angles A 0 , <i> 0 , and as indicated in figure 11. The angular orientation of the panel 



Figure 11.- Orientation of a panel relative to spacecraft center of mass. 

with respect to the spacecraft principal axes is given by 
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The torque exerted on the spacecraft by F is 


T = r x F 


(B4) 
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where 
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by application of equations (B3) and previous definitions. By using these results, the 
expansion of equation (B4) 


T x = yF z - zF y 
T y = “XF Z + zF x 
T z = xF y - yF x 


yields 
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after substituting equations (B2) and performing some lengthy reductions. Equations (B5) 
are the desired general expressions for the torque exerted on the spacecraft by a surface 
pressure acting on a panel having any orientation. 


Application to Hexagonal Spacecraft 

As all of the spacecraft panels were regular polygons having edges either parallel 
or normal to the spin axis z, the orientations of the x'-, y*-, and z'-axes for each panel 
can be chosen so that the angle A 0 is always zero. Equations (B5) then reduce to 

T x = -Fn(y sin ? 0 - z cos sin <£ 0 + y' sin cos 4> 0 - z' sin $ 0 ) 

+ F T |zfcos 3> 0 cos v - sin $ 0 sin £ 0 sin vj - cos sin v + cos (y' sin v - z' cos v) 


T y = f n( x sin ? 0 ' z cos ?0 cos $0 " y’ sin ?o sin ^o - z ' c o s %) 


x sin v - sin 


*o (y' 


+ jz^sin cos v + cos sin sin v} + cos 

— 

T z = -F n cos ? 0 (x sin 3> 0 - y cos *o " Y') ' F T x(cos <£ 0 cos v - sin <i> 0 sin C 0 sin v) 
+ y^sin cos v + cos $ 0 sin £ 0 sin vj + sin ? 0 ^y’ sin v - z' cos vj 


sin v - z’ cos v 


} (B6) 


For convenience in evaluating equations (B6), the side and end panels of the spacecraft 
were treated separately. 


Side panels .- As illustrated in figure 1, the panels comprising the sides of the 
hexagonal cylinder are parallel to the spacecraft spin axis so that the angle is zero 
for all six panels. With the spacecraft center of mass located as shown in figure 1, the 
origins of the x’~, y'-, and z'-axes for each of the six panels lie in the x-y plane at the 
constant distance iw/3 from the spin axis z such that 

Li 


r = 


/|w/3 cos 4> 0 \ 

i 

■|w\/3 sin 4> 0 | 

Li 

\° / 


Upon substituting these components of r, equations (B6) become 

T x = z'F n sin 4> 0 - F T j^w^3 sin sin v + cos 4> 0 (y' sin v - z’ cos 
Ty = -z’Fjq- cos + F-p|^w\/3 cos $ 0 sin v - sin <f> 0 ^y' sin v - z' cos v ) 
T z = y' F N “ | w /3F t cos v 
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End panels .- Referring again to figure 1, the ends of the hexagonal cylinder are 
normal to the spin axis z so that as defined in figure 11 is 90° and 270° for the 
solar cell and blank end panels, respectively. In this case the components of r are 


r = 


/ °\ 

I o ! 

1 

\ 2 / 


and equations (B6) reduce to 

T x = FN^Fy’ cos 4> 0 + z’ sin + 

T y = F]sf^y' sin 4? 0 - z' cos + 

T z = ?F^y' sin v - z' cos V s ) 

where the upper and lower of the double 
respectively. 


- ^ 
^■bF-p^tcos v cos 4> 0 - sin v sin 4> 0 ) 

ibF'j'^tcos v sin 4> 0 + sin v cos 4? 0 ) > 


(B7b) 


signs apply to the solar cell and blank end panels, 


Solar Radiation Pressure Torque 

In terms of the notation employed in figure 10 and equation (Bl), the solar radiation 
pressure force model of reference 8 is 


Fn = kj cos 77 + k2 cos 
Fp = k3 sin tj cos 77 



where 

k! = |p(l - s)SA 
k2 = (1 + sp)SA 
k 3 = (1 - sp)SA 

in which S is the intensity of the solar radiation pressure (see table Et) and 


(B8) 
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p fraction of incident radiation reflected by A 

1 - p fraction of incident radiation absorbed by A 

s fraction of p reflected specularly 

1 - s fraction of p reflected diffusely 

For the hexagonal configuration, the solar cell end panels were assumed to be per- 
fect absorbers of the incident radiation and all remaining panels, perfect specular 
reflectors. The factors p and s are therefore both unity for panels which are per- 
fect reflectors and both zero for those which are perfect absorbers. Hence, for each of 
the six side panels equations (B 7 a) reduce to 

T x = 2 z' SA cos^ 77 sin 

Ty = - 2 z’SA cos^rj cos (B 9 a) 

T z = 2 y'SA cos^ 77 v 

and for the end panels, equations (B 7 b) become 

T x = ^ SA sin 77 cos 77 cos 

Ty = | SA sin 77 cos 77 sin 

T z = 0 

where the resultants of y' and z' for both ends of the hexagon are zero for the space- 
craft orientation corresponding to the conditions listed in table HI, and 4 > 0 for each end 
become equivalent to arbitrary phase angles which may be set to zero. 

The expressions used to evaluate the angles 77 and v appearing in equations (B 9 ) 
were obtained by equating the direction cosines of the unit vector to the sun p in the 
same manner as for equations ( 4 ) and (A 3 ). The direction cosines of p relative to a 
spacecraft panel are given by equations (Bl) and by 
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1 cos L ' 


P = 


sin L cos e 


^ sin L sin ej 


(BIO) 


with respect to the geocentric inertial axes x^, y^, and Zj defined in figure 2 (also see 
ref. 11). By using equations (Bl) and (BIO) and the total transformation given by equa- 
tions (B3), (3), and (2), the desired relationships are 


/cos 77 \ 


1 cos L ' 

sin 77 cos v 

= c(x 0 ,$0,e 0 ) E H(«,r) 

sin L cos e 

ysin 77 sin vj 


^sin L sin e J 


In order to compute the torques given by equations (B9), the quantities y', z’, A, 
77 , and v must be determined. As a precise evaluation of this information for even a 
few spin cycles is not practical because of the large amount of computation required, the 
feasibility of making approximations was explored. Inspection of equations (Bll) shows 
that the main variations in y’, z T , A, 77 , and v are due to <p and ip, since the 
remaining angles are either constant or change very slowly. Hence, an obvious approach 
was to investigate the variations in the five quantities over a complete cycle of each <p 
and xj; with the rest of the angles in equations (Bll) held constant. 

Reference to figure 12, which illustrates the nominal initial spacecraft orientation 
relative to p for the assumed 45° sun- synchronous orbit, indicates the general manner 
in which the spacecraft panels experience solar radiation pressure. As this orientation 
will essentially be maintained by attitude control, the blank end panel will always be 
shielded from radiation pressure such that no torque can arise. The solar cell end will 
always be fully exposed to the sun, and the side panels will experience partial shielding 
by the solar cell panels in a periodic manner as the spacecraft spins. 

By using a partly computerized semigraphical technique for evaluating the six sets 
of y', z', and A for the side panels, the resultant torques given by equations (B9a) and 
(B9b) were found to be almost unaffected by small variations in p caused by the space- 
craft coning motion. As the resulting small fluctuations essentially averaged out, for 
coning angles within the 5° limit set by attitude control, 8 was set to zero with cp and 
xf/ replaced by their sum <J>. The error introduced by this approximation amounted to 
only a few percent for the hexagonal configuration and was judged to be negligible con- 
sidering the smallness and net effect of the torque itself. 
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y 



With the angles £ and L nominally constrained to differ by 45° and t bounded 
by attitude control limits, the only other effect to consider was the gradual change in the 
relative elevation of p from the initial 7.38° position shown in figure 12. This varia- 
tion, which annually increases and decreases the elevation of p by the obliquity of the 
ecliptic e, was found to cause a phase shift of the torque components accompanied by 
minor amplitude changes. The amount of phase shift proved to be equal to the projection 
of the difference between the angular separations of the spin axis z and p from the 
z^-axis on to a plane normal to the spin axis z. With this difference denoted by Ap, the 
resulting phase angle /3 is given by 


tan j3 = 


tan Ap 
cos (£ - L) 


(B12) 


The value of j3 for the initial Ap of 7.38° shown in figure 12 was about 10.4°. 

Substitution of $ for the sum of <p and \j/, where j8 is equivalent to sort of a 
composite 4> Q for the entire panel arrangement, permits reducing equations (Bll) to 
expressions for 77 and v which are sinusoids of the angle 4> - / 3 . Evaluation of equa- 
tions (B9) with these simplifications showed the resultant torque components to be mainly 
functions of <3? - j3 which appeared to be well approximated as 
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T x = -Rx sin ($ - /3) 
T y = -R y cos ($ - |3) \ 
T z = R z sin 6(4> - 0) 


(B13) 


where R x , Ry, and R z represent the resultant amplitudes of the SA-factors in equa- 
tions (B9a) and (B9b). Although y’, z', and A contained the expected harmonic at 6$ 
in all three members of equations (B9a), T x and Ty were dominated by equations (B9b) 
to the extent that only T z exhibited a distinct amplitude at 6#. 


Aerodynamic Pressure Torque 

The force model used to evaluate the aerodynamic pressure torque was that of ref- 
erence 9 and may be expressed as 


F N = kj cos 2 ri 


Fiji = k 2 sin rj cos 77 


(B14) 


where 


k l = 2(2 - cx N )q 0 A 
k 2 = cr T q Q A 


% = k p o v o 2 

in which A retains the same meaning defined in figure 10, and the surface reflection 
coefficients and ct-j- are both usually taken to be 0.8 as explained in reference 9. 

By proceeding in the same manner as before, equations (B7a) for the six side 
panels become 


T x = 0.80^. j 3 z' cos^ 77 sin <J ? 0 - sin tj cos t; jw /3 sin cos 1/ + 2 cos * 0 (y- sin v - z- cos z,)]J 
Ty = - 0 . 8 q o Aj 3 z' cos2tj cos 4> 0 - sin 77 cos 77 jw /3 cos $ 0 sin v - 2 sin <J > 0 (y' sin v - z' cos v ) | 


\ (B15a) 


T z = O-S^A^y' cos^ 7} - w /3 sin 77 cos rj cos v) 
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and equations (B7b) for the end panels reduce to 


T x = O-Sq^Ajs cos2 77 ^y' cos 4> 0 + z' sin <3> 0 ) + b sin 77 cos 

Ty = 0.8q o A|3 cos^ 77 (=Fy' sin 4> 0 - z' cos + b sin 77 cos 

T z = Tl.6q 0 A sin 77 cos 77 ^y' sin v - z' cos 


77 ^tcos y cos <£ 0 
77 ^tcos y sin <t> 0 


- sin y sin <I> 0 
+ sin v cos <i> 0 



(Bl5b) 


where the sign convention for equations (B7b) also applies to equations (Bl5b). 

With the direction of the incident aerodynamic pressure p assumed to coincide 
with that of the spacecraft velocity along the orbital path, the direction cosines of p 
relative to the x 0 -, y 0 -, and z 0 -axes defined in figure 2 are 


P = 


l°\ 


w 


(B16) 


Expressions for evaluating 77 and v in equations (B14) were then obtained from equa- 
tions (B16) and (Bl) by using the transformation given by equations (B3), (3), (2), and (1). 
The resulting relationships are 


cos 77 

sin 77 cos v 

= C(A 0 ,4> 0 ,? 0 ) E H(£,t) W(n,i,u)" 1 

l°\ 

1 

ysin 77 sin v J 


w 


(B17) 


As with equations (B9), the five quantities y*, z', A, 77 , and v must be known to 
calculate the aerodynamic pressure torques given by equations (B15). The variations in 
these quantities were evaluated in essentially the same manner as before, and similar 
approximations were made. However, this task was complicated in two ways by the 
spacecraft angle of attack. First, because the spacecraft will experience aerodynamic 
pressure at both positive and negative angles of attack, separate evaluations for each 
orientation mode were necessary. Secondly, the variation in torque amplitude with angle 
of attack had to be determined. 

Figure 13 depicts the spacecraft orientation relative to p for a positive angle of 
attack. For this case, the side panels will be partially shielded by the solar cell panels 
in a periodic manner similar to that for solar radiation pressure. The opposite situation 
occurs for negative values of the angle of attack a Q , as the solar cell panels then experi- 
ence periodic partial shielding by the side panels. 
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X 



Figure 13.- Spacecraft orientation relative to incident 
aerodynamic pressure. 

In accordance with the assumption in equation (B16) that p coincides with y Q of 
the spacecraft orbit, reference to figure 13 shows that a Q is determined by the scalar 
product of p with a unit vector along the spin axis z. By using the transformation 
between the x-, y-, and z-axes and x Q -, y Q -, and z Q -axes given by equations (1) to (3) as 

D = E(0,M/) m,T) W(ft,i,u)-1 

to obtain the components of z along the x Q -, y Q -, and z Q -axes the evaluation of this 
scalar product gives 

sin a Q = d 32 (B18) 

where d 02 is the second element in the third row of the D-matrix. 

The evaluations for the two orientation modes again showed the variations in y’, 
z', A, rj, and v to be well approximated with <p and xf/ replaced by their sum $ 
for coning angles up to 5°. For the case of positive a Q) the contributions from the side 
panels affected the total torque in much the same way as solar radiation pressure. 
Because of the particular choice for the panel dimensions, the torques arising from the 
blank end and side panels exactly canceled each other for negative values of q: 0 . The 
net torque was then entirely due to aerodynamic pressure on the back sides of the solar 
cell panels. The resultant torque components for each mode again were smooth sinu- 
soids of 4>. Because p rotates at the orbital frequency, an equivalent elevation change 
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is produced such that the angle u contributes a phase effect somewhat analogous to /3 
in equations (B13). 

The variation in torque amplitude with a Q was determined by evaluating equa- 
tions (B15) and (B17) over one cycle of $ for positive and negative values of a Q in 1° 
increments covering a range of 5°. Because of the approximation $ = $ + ip, with the 
coning angle set to zero, the values for a Q had to be obtained by incrementing t. This 
procedure showed the torque amplitudes to be essentially linear functions of a Q over 
the 5° range as expected. The resultant torques were surprisingly similar for both 
orientation modes, except that the z-components were of opposite sign. The respective 
magnitudes of T x and Ty for each mode differed by less than 2 percent, while those 
for T z were roughly within 5 percent. These results permitted replacing equa- 
tions (B15) with one set of expressions which apply to both modes. 

The resulting approximate equations for Tx> Ty, and T z were obtained by 
assuming £ = fl, as intended in table III, to permit writing the matrix product HW“lp 
in equation (B17) as 


1 

/»' 


j- sin u 

H(£,t) W(fl,i,u)-1 

1 

= 

cos u cos 


\°l 


^cos u sin cHq ^ 


where = i - r is the effective angle of attack for zero coning angles. Upon sub- 
stituting this result and approximating cos Q!q by unity, equations (B15) reduce to 

T x = A x Jsin <Xq j cos u cos ($ - u) > 

Ty = -Ay jsin a 0 | cos u sin ($ - u) \ (B19) 

T z = — Arj sin ot Q cos u sin 6 (# - u)^ 

where a Q ’ is now replaced by a Q which will represent the torque more accurately for 
situations when coning contributes an appreciable part of the effective angle of attack. 

The absolute value of sin a Q is necessary for the first two of equations (B19) since the 
signs of these components do not change with that of a Q as explained previously. 
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TABLE I.- SPACECRAFT CONFIGURATION AND TORQUE PARAMETERS 


Panel dimensions: 

Height, b, cm 101.60 

Width, w, cm 68.58 

Length, l, cm 111.76 

Moments of inertia and spacecraft weight: 

Transverse, I x = Iy = I, g-cm2 2.3885 X 10^ 

Axial, I z , g-cm2 2.6401 X 10^ 

Approximate ratio, Ig/I 1.2 

Approximate launch weight, kg 317.5 

Environmental torque coefficients: 

Residual magnetic moments, M x = My = M z , A(turn)-m2 0.07 

Eddy current, c, m^/ohm 1938.8 

Solar radiation pressure: 

Rx, dyne-cm 35.532 

Ry, dyne-cm 31.591 

R z , dyne-cm 0.678 

Aerodynamic pressure: 

A^, dyne-cm 338.87 

Ay, dyne-cm 382.60 

A;g, dyne-cm 24.99 
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TABLE II.- PHYSICAL PARAMETERS AND GENERAL CONSTANTS 
[From ref. 9, except as noted) 

Geodetic and geodynamical: 

Equatorial radius, r^, km 6378.2 

Mean radius, r a , km a 6371.2 

Universal gravitational constant, /x, km^/sec^ 398 601.5 

Axial rotation rate, co^, rad/sec 7.292115 x 10“® 

Obliquity of ecliptic, e, deg 23°27' 

Atmospheric and solar: 

Mean atmospheric density at 500 km, p Q , kg/m^ k 3.984 x 10“^^ 

Aerodynamic pressure at 500 km, q 0 = ^P 0 v o^» dynes/cm^ c 2.4134 x 10“^ 

Intensity of solar radiation pressure, S, dynes/cm2 4.5 x 10"^ 

a Reference 10. 
b Reference 12. 

c Calculated using p Q and v Q . 
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TABLE III. - INITIAL SPACECRAFT ORBIT AND STATE VARIABLES 


Elements for 500 -km. 45° sun-synchronous orbit: 

Inclination, i, deg 97.38 

Right ascension of ascending node, fi, deg 45.00 

Argument of latitude, u, deg 0.00 

Eccentricity, e 0 

Altitude, a, km 500.00 

Geocentric distance, r G = r^ + a, km 6878.2 

Mean motion, n 0 = rad/sec 0.001106798641 

Linear velocity along orbital path, v Q , km/sec 7.6164 

State variables: 

Angular momentum magnitude, h, dyne-cm-sec 27.9523 x 10^ 

Right ascension of K, deg 45.00 

Codeclination of h, t, deg 97.72 

Precession angle, 0, deg 0.00 

Coning angle, 6, deg 0.66 

Proper rotation angle, \p, deg 0.00 

Spin rate, u> z , rpm 3 

Attitude control: 

Maximum angle between orbit normal and spin axis, deg 5 
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